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Abstract 


>  /T/-  0  -J.  ^  ^  ^ 


-We-stu^jr  the  nonlinear  stability  of  discontinuous  steady  states  of  a  model  initial¬ 
boundary  value  problem  in  one  space  dimension  for  incompressible,  isothermal  shear  flow 
of  a  non-Newtonian  fluid  driven  by  a  constant  pressure  gradient.  The  non-Newtonian 
contribution  to  the  shear  stress  is  assumed  to  satisfy  a  simple  differential  constitutive  law. 
The  key  feature  is  a  non-monotone  relation  between  the  total  steady  shear  stress  and  shear 
strain-rate  that  results  in  steady  states  having,  in  general,  discontinuities  in  the  straiin  rate. 

..Yfo  showh^that  every  solution  tends  to  a  steady  state  as  t  and  we  identify  steady 


states  that  are  stable. 
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1.  Formulation  and  Discussion  of  Model  Problem 
We  study  the  quasilineax  system 


ut  =  Sx  , 
crt  +  <J  =  9{vx)  , 

on  [0, 1]  X  [0,  oo),  where 

S':=T  +  /x, 

with  /  a  fixed  positive  constant  throughout.  We  impose  the  botindary  conditions 


(1.1) 

(1.2) 


(1.3) 


S(0,t)  =  0  and  t;(l,t)  =  0, 


and  the  initial  conditions 


(1.4) 


accordingly, 


u(x,  0)  =  uo(2^)5  o'(x,0)  =  cro(x),  0<x<l; 


5(x,  0)  =  5o(x)  :=  <To(x)  +  uor(x)  +  fx  . 


(1.5) 


(1.6) 


The  function  g  ;  IR  — ♦  IR  is  assumed  to  be  smooth,  odd,  and  ^  >  0,  for  ^  0. 

The  initial-boimdary  value  problem  (1.1)-(1.5)  serves  as  a  simple  model  for  studying 
the  dynamic  behavior  of  a  non-Newtonian  fluid  in  Poiseuille  shear  flow  between  parallel 
plates  located  at  x  =  ±1;  the  flow  is  driven  by  a  constant  pressure  gradient  /.  In  this 
context,  u,  the  velocity  of  the  fluid  in  the  channel,  and  T,  the  shear  stress,  are  connected 
through  the  balance  of  linear  momentum  (1.1).  The  shear  stress  T  is  decomposed  into 
a  non-Newtonian  contribution  <7,  evolving  in  accordance  with  the  differential  constitutive 
law  (1.2),  and  a  viscous  contribution  Uj.  The  flow  is  assumed  to  be  symmetric  about  the 
centerline  of  the  channel.  Symmetry  and  compatibility  with  boundary  conditions  dictate 
the  following  restrictions  on  the  initial  data: 


uo(l)  =  0,  5o(0)  =  0,  and  <7-0(0)  =  0, 

which,  together  with  Eqs.  (1.2)-(1.4),  guarantee  that 


(1.7) 


I  COP'' 

<7(0,t)  =  Ux(0,f)  =  0 ,  (1.8) 

and  symmetry  is  preserved  for  all  time.  - 

The  system  (1.1)-(1.4)  admits  steady  state  solutions  (u(x),a(x))  satisfying  ^ - 

5  :=  (/(ux) -t- -f/x  =  0,  a  =  g(vr)  (1-9)  t 

■  _ 

on  the  interval  [0, 1].  In  case  the  function  w(^)  g(^)  -f  is  not  monotone,  there  may  be  - - 

multiple  values  of  Ui(x)  that  solve  Eq.  (1.9)  for  some  x's,  thus  leading  to  steady  velocity  - - - 

profiles  with  kinks  (see  Figs.  1  and  2).  Our  objective  is  to  study  the  stability  of  such  - - - 

ity  Codes 

_2_  and/or 

joist  I  Special 


□  D 


velocity  profiles;  we  also  study  well-posedness  and  the  convergence  of  solutions  of  (1.1)- 
(1.5)  to  steady  states  as  t  oo. 

Problem  (1.1)-(1.5)  captures  certain  key  features  of  a  class  of  viscoelastic  models  that 
have  been  proposed  to  explain  the  occurence  of  “spurt”  phenomena  in  non-Newtonian 
flows.  The  phenomenon  of  spurt  was  apparently  first  observed  by  Vinogradov  et  al.  [18] 
in  the  flow  of  highly  elastic  and  very  viscous  non-Newtonian  fluids  through  capilleiries 
or  slit  dies.  It  is  associated  with  a  sudden  increase  of  the  volumetric  flow  rate  occuring 
at  a  critical  stress  that  appears  to  be  independent  of  the  molectilar  weight.  It  has  been 
proposed  by  Himter  and  Slemrod  [5]  and  Malkus,  Nohel,  and  Plohr  [8],  [9],  and  [10]  that 
spurt  phenomena  may  be  explained  by  constitutive  laws  that  lead  to  a  nonmonotone 
relation  of  the  total  steady  shear  stress  versus  the  shear  strain-rate.  In  this  framework, 
the  increase  of  the  volumetric  flow  rate  corresponds  to  jumps  in  the  strain  rate  when  the 
driving  pressure  gradient  exceeds  a  critical  value. 


w(^)  X 


Fig.  1:  w  vs.  Fig.  2:  Velocity  profile  with  a  kink; 

w{-Vr{x))  =  fx. 

Interestingly,  various  three-dimensional  differential  constitutive  laws,  when  restricted 
to  one-dimensional  shearing  motions,  produ  ^  "onmonotone  steady  shear  stress  vs.  steady 
shear  strain-rate  relations.  The  simplest  such  model  leads  to  the  initial-boundary  value 
problem 

OCVt  ~  ^ X  ^^xx  f  ^ 

<Tt  -1-  <7  =  {Z  -|-  1)  Ur  >  (^•^^) 

Zt  Z  =  — <7Ur  • 

with  boundary  conditions  (1.4)  and  appropriate  initial  conditions;  the  parameters  a  and 
e  represent  the  ratio  of  Reynolds  number  to  Deborah  number  and  a  ratio  of  viscosities, 
respectively.  The  evolution  of  <7,  the  polymer  contribution  to  the  shear  stress,  and  of  Z, 
a  quantity  proportional  to  the  normal  stress  difference,  are  governed  by  the  second  and 
third  equations  in  (1.10),  which  are  a  restriction  of  differential  constitutive  laws  due  to 
Oldroyd  [12]  and  Johnson  and  Segalman  [6]  (see  [16]  or  [8]  for  a  derivation).  The  steady 


state  solutions  (u(x),a(x),  Z(x))  satisfy 

(1.11) 

and  the  steady  strain-rate  is  determined  by  solving  the  equation 

T(Vj:)  ■•=  —7^2"  +£Ux  =  -/g  •  (1-12) 

1  -f-  u* 

For  £  <  |,  T{i)  is  not  monotone;  its  graph  has  the  shape  shown  in  Fig.  1.  Malkus,  Nohel 
and  Plohr  [S]  performed  numerical  calculations  for  solutions  of  (1.10)  that  reproduce  the 
experimental  results  on  spurt  of  Vinogradov  et  aJ.  [18]  qualitatively  and  quantitatively. 
Moreover,  motivated  by  the  relative  scales  of  the  parameters  a  and  e  in  the  experiments, 
they  study  an  associated  system  of  ordinary  differential  equations  that  is  obtained  through 
a  “quasidynamic  approximation”  (setting  or  =  0)  in  (1.10).  Their  analysis  of  the  approxi¬ 
mating  system  in  [9]  and  [10]  predicts  spurt  as  well  as  latency,  shape  memory  and  hysteresis 
under  cyclic  loading. 

Here,  we  test  the  hypothesis  that  lack  of  monotonicity  of  the  steady  shear  stress 
function  causes  non-Newtonian  response,  at  the  level  of  the  paradigm  (1,1)-(1.5).  We 
remark  that  this  problem  can  be  regarded  as  a  formed  approximation  to  (1.10)  arising 
by  fixing  Z  in  (I.IO2)  at  its  steady  state  value  in  (1.11).  For  simplicity,  the  function 
=  ^(0  +03  assumed  to  have  one  single  loop.  The  loop  forms  as  the  combined 
outcome  of  the  non-Newtonian  contribution  to  the  steady  shear  stress  associated  with  g{^) 
and  the  Newtonian  contribution  associated  with  The  hypotheses  on  g(^)  imply  that 
to(if)  is  odd  and  tu(^)  ^  0  for  ^  7^  0.  The  graph  of  a  representative  w[^)  for  ^  >  0  is  shown 
in  Fig.  1;  in  the  figure,  m  and  M  stand  for  the  levels  of  the  bottom  and  top  of  the  loop, 
respectively.  Our  smalysis  can  be  routinedly  generalized  to  cover  the  case  when  iy(0 
a  finite  number  of  loops,  but  we  do  not  pursue  this  direction  here. 

Steady  state  solutions  for  the  representative  w{^)  are  constructed  as  follows:  First, 
solve  k;(u(x))  =  fx  for  each  given  x  €  [0,1],  where  u(x)  =  —Vx{x).  This  equation  admits 
a  unique  solution  for  0  <  /x  <  m  or  /x  >  M,  and  three  solutions  for  m  <  fx  <  M.  Let 
V(x),  0  <  X  <  1,  be  such  a  solution.  Setting 

v{x)  =  J^u{y)dy,  a{x)  =  g{-xl{x))  (1.13) 

(v(x),a(x))  satisfy  (1.9)  and  (1.4)  and  give  rise  to  a  steady  state.  Clearly,  if  /  <  m  there 
is  a  unique  smooth  steady  state,  li  m  <  f  <  M  there  is  a  unique  smooth  steady  velocity 
profile  and  a  multitude  of  profiles  with  kinks,  finally,  if  /  >  M  all  steady  velocity  profiles 
have  kinks.  An  example  of  a  velocity  profile  with  a  kink  is  shown  in  Fig.  2. 

We  are  interested  in  studying  the  dynamics  of  solutions  to  (1.1)-(1.5),  and  in  partic¬ 
ular,  to  determine  which  of  the  steady  states  are  stable.  We  establish  that  every  solution 
of  (1.1)-(1.5)  converges  as  <  00  to  a  steady  state  solution  (u(x),a(x)).  It  would  be 

interesting  to  identify  the  region  of  attraction  of  the  stable  steady  solutions,  however,  this 


v. 


a  = 


l  +  vl' 


Z  +  l 


1  -}-  u 


5t2  ’ 
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seems  a  very  complicated  task.  Instead,  we  show  that  steady  states  (v(x),a(x)),  with  a 
finite  number  of  jumps  in  Vx{x)  or  a(x),  and  such  that 

wX~Vx{x))  >  Co  >  0,  a.e.  x  €  [0,1]  (1-14) 

are  “asymptotically  stable”  in  a  sense  that  is  made  precise  in  Theorem  5.1.  Note  that 
relation  (1.14)  excludes  solutions  that  take  values  in  the  monotone  nonincreasing  part 
of  the  graph  of  w(^).  The  stable  solutions  are  local  minimizers  of  an  associated  energy 
functional. 

The  above  results  are  close,  in  spirit  and  in  technique,  to  the  analysis  of  Andrews 
and  Ball  [1]  and,  especially,  that  of  Pego  [14]  for  phase  transitions  in  one-dimensional 
viscoelastic  materials  of  rate  type.  Current  work  of  Novick-Cohen  and  Pego  [11]  on  spinodal 
decomposition  involves  similar  ideas. 

The  paper  is  organized  as  follows:  In  Sec.  2, -we  obtain  basic  a-priori  estimates  and 
determine  invariant  regions  for  an  associated  ordinary  differential  equation.  In  Sec.  3, 
we  discuss  the  well-posedness  and  regularity  of  solutions  of  a  system  incorporating  (1.1)- 
(1.5),  a  consequence  of  which  are  existence  and  regularity  results  for  (1.10)  and  other 
popular  models  for  motions  of  non-Newtonian  fluids  arising  from  differential  constitutive 
laws.  In  Sec.  4,  we  study  the  convergence  of  classical  solutions  of  (1.1)-(1.5)  to  steady 
states.  Finally,  in  Sec.  5,  we  state  and  prove  the  main  result  on  stability  of  steady  states 
by  making  use  of  the  analysis  in  Secs.  2  and  4. 


2.  A  Priori  Estimates  and  Invariant  Regions 

Problem  (1.1)-(1.5)  is  endowed  with  certain  identities  that  play  a  crucial  role  in  the 
analysis.  As  a  consequence  of  (1.1)  and  (1.4),  smooth  solutions  of  the  system  (1.1)-(1.5) 


satisfy: 

1/2—  /  v^dx+  1  Svxdx  =  0  , 
d't  Jo  Jo 

(2.1) 

and 

d 

1/2—  /  Vfdx  +  /  dx  +  /  (St  +  S)vxtdx  —  0 . 

dt  Jq  Jq  Jq 

(2.2) 

Using  (1.3)  and  (1.2),  Eq.  (2.2)  leads  to  the  energy  identity 


^{1/2  J  v1dx  +  J  +  X  fvj;]dx}  +  J  [u^ -f- =  0  .  (2.3) 


The  function  ^  ^ 

W(i)  :=  I  w(0<i<:  =  ine  +  /  g(OdC,  (2.4) 

JQ  Jq 

plays  the  role  of  a  stored  energy  function,  and  is  not  convex.  This  fact  is  the  main  obstacle 
in  the  analysis  of  stability.  Since  ig{i)  >  0,  it  follows  that  g{Od(  >  0  for  and 

IF(^)  satisfies  the  lower  bound 

>^(0  +  fx^  >  -  f  ,  eelR,  0  <  I  <  1 .  (2.5) 
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Integrating  (2.3),  one  has  the  identity 


1/2  f  v1{x,t)dx+  f  {W{vx{x,t))  +  fxvx{x,t))dx  +  f  f  +  vli)dxds  = 

Jo  Jo  Jr  Jo 

1/2  I  Vt{x,T)dx+  j  (W(vx{x,t)  +  fxVx{x,T))dx,  0<T<t. 

Jo  Jo 


(2.6) 


As  a  first  consequence,  the  energy  identity  (2.6)  yields  an  a  priori  bound  for  S  as 
follows.  First  observe  that  (2.5)  and  (2.6)  with  r  =  0  imply  the  estimate 


f  v'}(x,t)dx  +  1/2  /  vl(x,t)dx  +  2  /  f  {v'f  +  vl^)dxds 
Jo  Jo  Jo  Jo 

<  2p  +  [  SQ^(x)dx  +  2  f  [W(t;oi(x))  +  fxvQx{x)]dx  <  C , 
Jo  Jo 


(2.7) 


where  C  is  a  constant  depending  only  on  the  data.  From  (1.1)  and  (1.4)  we  obtain  the 
inequality 


5^(x,t)  =  [/  U4(^,<)d^]^  <  /  v^(x,t)dx, 
Jo  Jo 


which,  together  with  (2.7),  implies 


|5(x,t)|<C  0<x<l,  0<i<oo, 


(2.8) 


(2.9) 


Control  of  S  enables  us  to  take  advantage  of  the  special  structure  of  Eq.  (1.2).  It  is 
convenient  to  introduce  the  quantity 


s:=  <7  +  f  X. 

Then,  Eqs.  (1-2),  (1.3),  and  (2.10)  readily  imply 

St  +  s  +  g{s  —  S)  =  f  X. 


(2.10) 


(2.11) 


For  a  fixed  x,  Eq.  (2.11)  may  be  viewed  an  ODE  with  the  forcing  term  5(x,  •).  Also, 
observe  that  for  a  steady  state  (a,Vx),  one  has  5  =  0,  and  s  =  —Vx  is  an  equilibrium 
solution  of  (2.11)  (with  5  =  0).  If  5  =  0  in  (2.11),  it  is  evident  that  the  ODE  admits 
positively  invariant  intervals  for  each  fixed  x.  We  claim  that  this  property  is  preserved  in 
the  presence  of  a  priori  control  of  5.  Such  control  is  provided  by  bounding  the  L^-norm 
of  Vf  in  (2.8)  as  in  (2.9);  more  delicate  bounds  are  derived  in  Sec.  5. 

To  fix  ideas,  assume  that  for  0  <  t  <  to 


|5(x,t)|<p  ,  0<x<l 


(2.12) 


for  some  p  >  0.  For  x  fixed  in  [0, 1],  we  use  the  notation  5(t)  :=  5(x,t)  and  conveniently 
rewrite  (2.11)  as 

5,  +  w(3  -  S(i))  =  fx-  S{t) .  (2.13) 
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We  state  the  following  lemma  on  invariant  intervals;  its  proof  is  obvious. 

Lemma  2.1.  Let  S  satisfy  the  uniform  bound  (2.12)  for  Q  <t  <to.  For  x  fixed,  0  <  x  <  1, 
assume  there  exist  s_,  3+  such  that  s_  <  5+  and 

w{s^  —  X)<fx  —  X  ,  1A|  < /?  (2-14) 

w(s+  —  X)>fx-X  ,  |A|  <  p  (2.15) 

Then  the  compact  interval  [s_,  s+]  is  positively  invariant  for  the  ODE  (2.13)  on  the  time 
interval  0  <  t  <  <o- 

With  the  goal  of  generating  invariant  intervals,  we  study  the  solution  sets  of  the 
inequalities  (2.14)  and  (2.15)  as  fvmctions  of  p  and  x.  First  observe  that,  since  ^  lim  w(^)  = 

±co,  given  any  x  and  p  one  Ccin  determine  sq+  large,  positive  and  Sq-  large,  negative  such 
that  if  s_  <  So-  aJ^d  s+  >  so+,  then  s_  and  s^.  satisfy  (2.14)  and  (2.15),  respectively. 


W  (^) 


The  possibility  for  more  discriminating  choices  of  invariant  intervals  arises  if  we  restrict 
attention  to  small  values  of  p.  For  a  function  ty(^)  with  a  single  loop,  the  most  interesting 
case  arises  when  f  x  —  p,  f  x  and  f  x  -{■  p  each  intersects  the  graph  of  it^(0  three  distinct 
points.  Referring  to  Fig.  3,  the  abscissae  of  the  points  of  intersection  are  denoted  by 
(q_,;5_,7_),  (q;o,/3o,7o)  and  (a+,/3+,7+),  respectively.  First  examine  (2.15).  For  0  < 
X  <  p,  any  s+  €  {oiq  +  p,f3o)  U  (70  +  p,oo)  satisfies  (2.15),  while,  for  -p  <  A  <  0,  any 
S4.  €  -  p)  U  (74., 00)  satisfies  (2.15)  (here,  by  convention,  any  interval  that  is 

not  well  defined  will  stand  for  the  empty  set).  Thus,  the  solution  set  of  (2.15)  contains 
an  unbounded  component  Ju  ■=  (70  +  P>oo)  f"*  (7+>o°)>  well  as  a  bounded  component 
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Jb  '•=  (ao+p,^o)n(a+,;5+-p).  Correspondingly,  the  solution  set  of  (2.14)  dso  contains  two 
components  :=  (-oo,a_)n(-oo,ao -p)  and  h  ■=  (/?o,7o  -p)n(/?_  +p,7-).  Choosing 
s-  €  /«  and  s+  €  Ju  we  construct  invariant  intervals  [s_,s+]  of  the  type  discussed  in 
the  previous  paragraph.  More  interesting  choices  occur  if  one  chooses  p  so  small  that 
and/or  Jb  are  nonempty;  one  then  takes  s_  €  /&  and  s+  £  or  s_  G  and  s+  G  Jb  (see 
Fig.  3). 

The  above  discussion  identifies  invariant  intervals  for  (2.13)  when  x  is  fixed.  Of 
particular  interest  to  us  is  to  study  the  persistence  and  dependence  on  x  of  the  invariant 
intervals  as  x  varies  in  [0,1].  In  this  context  p  is  fixed  and  (2.12)  is  assumed  to  hold. 
lu,  Ju,  Jb  and  Jb  are  defined  as  above  and  depend  on  x  as  well  as  p. 

If  we  choose  s_  G  and  G  P|  (which  are  nonempty)  then  [s_,s+]  is 

*e(o,i]  i€{o,i] 

an  invariant  interval  for  (2.13)  valid  for  all  x  G  [0,  Ij.  As  an  application  of  this  remark  and 
by  viitue  of  the  bound  (2.9),  we  conclude  that 

|s(x,<)l<C,  0<x<l,t>0  (2.16) 

which,  in  turn,  using  (1.3)  and  (2.10),  implies 

|u;c(x,  t)|  <  C ,  0  <  X  <  1,  f  >  0  .  (2.17) 

By  contrast,  the  discriminating  invariant  intervals  degenerate  as  we  approach  the  top 
or  bottom  of  the  loop.  However,  if  p  is  small  enough  and  we  avoid  the  top  or  bottom  of 
the  loop,  the  invariant  intervals  may  be  chosen  of  uniform  length.  In  Sec.  5  we  use  a  result 
along  these  lines  which  we  formulate  below. 

T-et  7(1)  be  a.  piecewise  smooth  solution  of 

w(s(x))  =  fx  (2.18) 

defined  on  [0, 1]  and  admitting  jump  discontinuities  at  a  finite  number  of  points  Xi, ...,  Xn 
in  [0,  Ij.  In  addition,  suppose  that  s{x)  takes  values  in  the  monotone  increasing  parts  of 
the  curve  w{Q  and  that  it  avoids  jumping  at  the  top  or  bottom  of  the  loop,  i.e., 

w'{J(x))  >  cq  >  0,  X  G  (0, 1]  \  {xi,...,x„},  (2.19) 

for  some  constant  cq.  We  look  for  invariant  intervzds  for  (2.13)  that  are  of  uniform  length, 
centered  around  ^(x)  at  each  x.  First,  consider  the  case  p  =  0  in  (2.12),  i.e.  5  =  0, 
in  (2.13).  For  each  fixed  x,  let  i?o(3;)  be  the  distance  between  s(x)  and  the  nearest  solution 
of  the  equation  w(^)  =  fx  different  from  ^(x).  Set  i?o  :=  info<i<i  7?o(3;)  and  observe 
that  by  (2.19),  Rq  >  0.  Then  for  any  r  <  Rq,  [s(x)  —  r,  s(x)  +  r]  is  an  invariant  interval 
for  (2.13)  with  5  =  0.  For  p  >  0,  a  review  of  the  above  constructions  yields: 

Corollary  2.2.  Given  s(x)  satisfying  (2.18)  and  (2.19),  there  are  positive  constants  po 
and  Rq,  such  that  given  any  p  <  po  there  are  positive  numbers  rj[(p)  and  r2(p)  =  Rq  —  rj(p) 
with  the  property:  If  (2.12)  holds  on  [0,fo],  then  for  any  r  G  (ri(p),  r2(p)),  the  interval 
[s(x)  —  r,s(x)  +  r]  is  positively  invariant  for  the  ode  (2.13)  on  [0,<o]  for  any  x  G  [0,1]  \ 
{xi,...,x„}.  The  function  r\{p)  is  an  increasing  function  of  p,  which  is  of  the  order  of  p, 
as  p  — V  0. 
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3.  Existence  and  Regularity  Theory 

In  this  section,  we  study  existence  and  regularity  questions  for  a  coupled  system 
of  partial  differential  equations  that  describes  one-dimensional  shearing  motions  of  non- 
Newtonian  fluids.  This  system  is  appropriate  for  non-Newtonian  fluid  models  in  which  the 
total  stress  is  decomposed  into  a  Newtonian  contribution  2md  a  finite  number  of  relaxation 
stress  components,  viewed  as  internal  variables  and  evolving  in  accordance  with  differential 
constitutive  laws.  Some  specific  models  include  (1.10)  and  others  that  axe  discussed  in  the 
Appendix. 

Let  u  be  a  scalar  and  U  =  {Ui,  be  a  N-vector  and  consider  the  system  of 

partial  differential  equations 

Vt  —  Sx  ) 

Ut=H%U,Vx,x),  (3.1) 

S  =  Vx+G-(U,x)  , 

on  [0, 1]  X  [0,oo).  In  the  context  of  viscoelasticity,  u  is  a  scaJar  velocity  and  U  represents 
a  vector  of  relaxation  stresses  viewed  as  internal  variables.  For  definiteness,  we  impose 
boundary  conditions  appropriate  for  symmetric  channel  flow  (some  other  boundary  condi¬ 
tions  may  be  treated  in  a  similar  fashion); 

5(0,0  =0,  v(l,0  =  0,  '  (3-2) 

and  initial  conditions  of  the  form 

v(x,0)=v,(x),  U(x,0)=U.{x),  I  6  [0,1).  (3.3) 

We  assume  that  the  scalar  function  G*  6  and  the  vector  function  H*  €  with  respect 
to  their  arguments  for  U  €  Vx  €  R,  and  x  €  [0, 1].  We  indicate  in  the  text  whenever 
different  smoothness  hypotheses  are  imposed  on  the  functions  G*  and  H*. 

Two  parallel  existence  and  regularity  theories  will  be  pursued.  The  first  is  appropriate 
for  smooth  initial  data  and  yields  classical  solutions.  The  second  allows  discontinuities  in 
Vor(2^)  3^d  in  Uq{x)  and  yields  strong  solutions.  The  latter  result  allows  prescription  of 
discontinuous  initicd  data  of  the  same  type  as  the  discontinuous  steady  states  studied  in  this 
paper.  Finally,  we  develop  criteria  for  global  existence  which  apply  to  the  models  at  hand, 
namely,  to  systems  (1.1)-(1.5),  (1.10),  as  well  as  to  some  more  generaJ  models  discussed 
in  the  Appendix.  We  note  that  existence  of  classical  solutions  to  the  system  (3.1)-(3.3) 
may  also  be  obtained  by  using  an  approach  based  on  the  Leray  -  Schauder  fixed  point 
theorem  (cf.  Tzavaras  [17]  for  existence  results  on  a  related  system).  Other  existence 
results  for  (1.10)  were  obtained  by  Guillope  and  Saut  [2],  md  small-data  results  for  some 
other  models  in  more  than  one  space  dimension  in  [3]. 

Our  first  result  establishes  local-time  existence  of  classical  solutions  of  (3.1)-(3.3). 
Theorem  3.1  (Local  classical  solutions):  Assume  that  Uo(l)  =  0,  and  that  with  So{x)  = 
+  G*{Uo{x),x),  we  have  So  €  fr*[0,l]  for  some  s  >  3/2,  with  5<,(0)  =  0,  vo(l)  = 
5„(1)  =  0,  and  that  Uo  €  C'[0, 1]^,  where  denotes  the  usual  interpolation  space  [7]. 
Then  for  some  T  >  0  a  unique  solution  of  (3.1)-(3.3)  exists  with  v  €  G([0,  T),G^),  vt  € 
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C([0,T),C°),  U  €  T],  (C^)^).  IfT  is  the  maximal  time  of  existence,  butT  is  finite, 

then 

limsup[  sup  (|ui(x,<)|  +  [{/(Xji)!)]  =  oo  .  (3.4) 

t— T-  *e[o,i] 

Proof.  It  is  convenient  to  study  the  problem  with  v  replaced  by  S.  After  computing  St 
from  (3.1)  and  substituting  for  u*  =  (5  —  G*),  we  formedly  obtain  equations  of  the  form 

St=S,r+G(U,S,x),  Ut  =  H{U,S,x)  (3.5) 

where  G  and  H  are  functions.  S  should  satisfy  the  boundary  conditions 

5(0,t)  =  0,  5,(l,t)  =  0.  (3.6) 

We  study  existence  for  the  degenerate  semilinear  reaction-diffusion  system  in  (3.5) 
using  the  standard  approach  of  Henry  [4].  We  consider  (3.5)  ais  an  abstract  parabolic 
equation  zt  +  Az  =  f{z)  on  a  Banach  space  X.  Here,  z  =  (5,  U),  the  operator  A  is  defined 
by  A{S,U)  =  (-5xr,0)  and  the  Banach  space  X  =  1]  x  C^[0, 1]^.  The  domain  of 

the  Laplacian  —AS  =  —Sxz  in  X*  is 

D(-A)  =  {5|  SeH'^[0,l]  with  5(0)  =  0,  Sr(l)  =  0}  , 

so  the  domain  oi  A  is  D{- A)  x  .  The  Laplacian  is  a  sectorial  operator  in  L^,  so 

A  is  a  sectorial  operator  on  X,  and  generates  an  analytic  semigroup.  For  0  <  a  <  1,  the 
domaiin  of  the  fractional  power  A°  is  X°‘  =  £)((— A)“)  x  (C^)^ .  For  o  >  3/4  this  domain 
is  continuously  embedded  into  x  (C^)^,  where 

C,'  =  {5€C'[0,1]|  S(0)  =  0,  5,(1)  =  0}. 

The  hypotheses  on  G*  and  imply  that  the  map  {S,U)  — *•  {G{U,S,x),H{U,S,x))  is 
as  a  map  from  X^  to  X.  The  hypotheses  on  the  initial  data  imply  that  initially 
2(0)  e  X°  =  D{A°)  for  some  a  >  3/4.  The  proof  of  Theorem  3.3.3  of  Henry  [4]  applies, 
yielding  existence  of  a  unique  solution  2(t)  of  the  abstract  equation  satisfying,  for  this  a, 
and  for  all  7  <  1, 


2  G  C([0,T),X“)nC'((0,T),X^)nC((0,r),I?(A)).  (3.7) 

(See  Lemma  3.2.1  and  Theorem  3.5.2  in  [4]  concerning  differentiability  of  the  solution.  For 
an  alternative  source,  see  Pazy  [13].)  Taking  components,  we  find  that 

5  6  C([0,T),Q')nC‘((0,T),Cj')nC((0,T),i?(-A)),  i7e  Cn[0,r),C'); 

the  latter  follows  by  taking  t  to  zero  in  (3.5).  Now,  recover  v  from 

t>(i,i)  =  -  f  [S(y,t)-G-iU{v,t),y)ldy. 
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It  follows  that 


t>  e  C(10, T), C=)  n  C‘((0,T), C’),  V,  =  s,e  C((o, t), C)  , 

the  system  (3.1)  is  satisfied  for  <  >  0,  and  v{x,0)  =  vo{x). 

To  prove  the  assertion  regarding  blow-up  if  global  existence  fails,  we  apply  Henry’s 
Theorem  3.3.4  to  conclude  that  if  the  maximal  T  is  finite,  then  the  —  norm  of  the 
solution  (5,  U)  of  (3.5)  must  blow-up  as  t  +  T— .  But  then  it  follows  that  the  sup  norm  of 
(5,  U)  must  blow  up,  because  otherwise,  the  nonlinear  terms  f{z)  in  the  abstract  equation 
satisfy  a  linear  growth  estimate  ||/(-2)1|a'  ^  -^^(1  +  Ikllx^),  a^id  global  existence  follows 
from  Corolla-y  3.3.5  in  [4].  Now  the  assertion  in  (3.4)  follows  from  the  blow-up  of  the  sup 
norm  of  (5,  U).  | 

Strong  “semi-classical”  solutions  with  discontinuous  stress  components  may  be  ob¬ 
tained  by  a  different  choice  of  function  spaces  in  the  proof  above.  The  following  local 
existence  theorem  yields  solutions  of  (3.1)  in  which  U  and  Vx  may  be  discontinuous  in 
X,  but  Sx  and  Vt  are  continuous,  and  17  is  as  a  function  of  t  for  every  x.  Thus  all 
derivatives  in  (3.1)  may  be  interpreted  in  a  classical  sense  as  long  as  the  equation  is  kept 
in  conservation  form.  A  result  of  this  type  was  obtained  by  Pego  in  [14];  we  follow  a  similar 
line  of  argument.  In  what  follows  W^’°°  is  the  space  of  Lipschitz  continuous  functions  on 
[0, 1],  and  is  the  space  of  functions  with  Lipschitz  continuous  first  derivatives. 

Theorem  3.2  (Local  strong  solutions):  Assume  that  dG* /dU  and  H*  are  locally  Lipschitz 
functions.  Assume  that  Uo(l)  =  0)  So{^)  =  Vox{^)  +  G*(Uo(x),x),  we  have 

So  e  H^[0,1]  with  5o(0)  =  0,  and  that  Uq  €  L°°[0, 1]^.  Then  for  some  T  >  0  a  unique 
solution  of  (3.5)-(3.6)  exists  with 

S  eC{[0,T),H^)nC({0,T),W^’‘^)nC\{0,T),H^)  for  all  s<2, 

U  €C\[0,T),{L°^f). 

Moreover,  there  exists  a  unique  solution  of  (3.1)-(3.3)  with 

V  €  C([0,  T),  n  C'((0,  T),  . 

If  T  is  maximal  hut  finite,  then  (3-4)  holds.  Also,  given  a  hounded  representative  Uo-  of 
the  equivalence  class  Uo,  there  exists  T,  >  0  and  a  unique  hounded  measurable  function 
U»{x,t),  representing  U[-,t')  for  each  t,  such  that  t  — >•  C/,(x,t)  is  on  [0,T,)  for  each  x, 
and,  identifying  the  equivalence  class  S  with  its  unique  continuous  representative,  (S',  17,) 
is  a  strong  solution  of  (3.5). 

Proof  of  Theorem  3.2.  The  main  part  of  the  proof  is  a  slight  modification  of  the  proof  of 
Theorem  3.1.  Now,  the  Banach  space  X  =  L^[0, 1]  x  L°°[0, 1]^.  One  may  show  that  with 
a  =  1/2,  D{{-A)°‘)  =  n  {5  I  5(0)  =  0},  and  this  space  is  continuously  embedded  in 
C[0,1].  One  obtains  a  solution  z(<)  to  the  abstract  equation  with  the  properties  in  (3.7) 
for  a  =  1/2,  with  X°'  =  i7((-A)“)  x  L°®[0, 1]^.  The  regularity  properties  for  the  solu¬ 
tion  follow  by  taking  components  and  using  the  equations  satisfied  by  Ut  and  Sxx-  The 
assertions  regarding  blow-up  follow  eis  before. 


-11- 


To  establish  the  existence  of  the  ‘classical’  solution  Z7,,  we  argue  as  follows.  Identify 
the  equivalence  class  S  with  its  unique  continuous  representative,  and  for  each  fixed  x, 
let  U^{x,t)  be  the  unique  solution  of  the  ODE  dUt,/dt  =  H{U^,S,x)  with  iniliaJ  value 
Uomix).  Solutions  exist  on  some  interval  [0,T.),  T*  <  T,  with  T,  indenpendent  of  x,  but 
the  maximal  T*  may  in  principle  depend  on  Uo*  and  on  how  H  depends  on  x.  Considering 
Ut  =  H{U,S,x)  as  an  ODE  on  it  follows  from  uniqueness  that  Um  agrees  with  U 

for  almost  all  x,  for  all  f  6  [0,T*).  | 

In  the  rest  of  this  paper,  U  will  usually  stand  for  a  ‘classical’  solution  17,.  We  remark 
that  CIS  a  simple  corollary  of  the  continuous  dependence  of  solutions  of  dU^,/dt  =  H{U»,  5,  r) 
on  initial  data  (or  final  data),  it  follows  that  the  solution  17,  is  continuous  at  (io,to)  if 
and  only  if  Uo*  is  continuous  at  Xo-  In  particular,  discontinuities  in  the  initial  data  for 
U  are  preserved.  In  general,  the  spatial  smoothness  of  the  solution  of  (3.1)-(3.3)  may  be 
limited.  However,  so  long  as  G*  and  H*  are  smooth,  the  solution  is  smooth  as  a  function 
of  time:  The  following  proposition  is  a  straightforward  application  of  Corollary  3.4.6  in  [4]. 
The  hypotheses  guarantee  that  G  and  H  are  functions  of  (5, 17),  hence  the  nonlinearity 
f[z)  in  the  abstract  equation  is  C’’  as  a  map  from  X°  to  X.  Henry’s  Corollary  3.4.6  in  [4] 
asserts  that  the  mapping  t  — +  z{t)  is  C’’  as  a  map  on  (0,  T)  with  values  in  X°‘. 
Proposition  3.3  (Temporal  regularity):  Assume  that  partial  derivatives  ofG*  with  respect 
to  U  up  to  order  k  +  1  are  continuous,  and  partial  derivatives  of  H*  with  respect  to  U  and 
Vx  up  to  order  k  are  continuous. 

a)  Suppose  r  =  fc  —  1  >  1.  Under  the  hypotheses  of  Theorem  3.1,  we  have 

S  €  ^((0, T), Ci),ue  C'((0,T), (C'f),vs  ^((0, T), C^) 

b)  Suppose  r  =  fc  >  1.  Under  the  hypotheses  of  Theorem  3.2,  we  have 

S  €  C"((0,  T),H^),Ue  ^((0,  T),  ,  v  e  C’'((0,  T),  W'’~) . 

We  conclude  this  section  by  discussing  criteria  that  imply  the  global  existence  of 
solutions  of  (3.1)-(3.3).  Under  the  assumptions  of  Theorem  3.1  (or  Theorem  3.2),  the 
initial-boundary  value  problem  (3.1)-(3.3)  has  a  unique  solution  (t;(r,t),  U{x,t))  on  [0, 1]  x 
[0,T)  for  some  (maximal)  T  >  0.  In  view  of  (3.4)  ,  if  Vx{x,t)  and  U{x,t)  are  a-priori 
bounded  on  any  interval  [0, 1]  x  [0,  to]>  with  to  >  0,  then  T  =  oo  and  the  solution  is  defined 
for  all  times.  In  light  of  (2.16)  and  (2.17)  this  is  the  case  for  the  problem  (1.1)-(1.5). 
Corollary  3.4.  Under  the  hypotheses  stated  in  the  Introduction,  the  initial-houndary  value 
problem  (1.1)-(1.5)  has  a  unique  solution  existing  globally  in  time. 

For  a  class  of  problems  describing  viscoelastic  fluid  flows,  including  (1-11)  and  others 
considered  in  the  Appendix,  the  function  H*  in  (3.1)  has  linear  growth  in  u*.  For  certain 
of  thes^  -"odels  it  can  be  shown  that  U  satisfies  an  a  prion  L°°  estimate  on  [0, 1]  x  [0,<o] 
for  a  .v  ,  >  0  (cf.  Appendix),  and  this  suffices  to  guarantee  that  solutions  are  globally 
define  V... 

Theorei*^  (Global  Existence).  Let  the  assumptions  of  Theorem  3.1  (or  Theorem  3.2) 
he  s  tisfied  and  let 

H-{U,Vx,x)  =  Hi{U,x)vx  +  H2{U,x).  (3.8) 


-12- 


Lei  (v(x,i),  U{x,t))  be  a  solution  of  (S.1)-(3.S)  as  in  Theorem  S.l  (or  Theorem  3.2) 
defined  on  [0,1]  x  [o,T),  for  some  T  >  0.  If  T  is  the  maximal  time  of  existence,  hut  T  is 
finite,  then 

linisup[  sup  |C/'(x,<)|]  =  c» .  (3-9) 

t— T-  re[o,i] 

Proof  of  Theorem  3.5.  Let  (u(x,<),  U{x,t))  be  a  solution  of  (3.1)-(3.3)  on  [0,1]  x  [0,T), 
■with  T  maximal,  and  assume  that  (3.9)  is  violated;  that  is,  there  is  a  positive  constant  Ki 
such  that 

\Uix,t)\<Ki  for  (x,t)6[0,l]x[0,T).  (3.10) 

In  view  of  (3.1)  and  (3.4),  it  suffices  to  show  that 

\S{x,t)\<K2  for  (x,<)g[0,1]x[0,T),  (3.11) 

where  K2  is  some  finite  constant. 

Indeed,  by  (3.1)  and  (3.8),  S(x,t)  satisfies  a  linear  parabolic  equation  of  the  form 

S,  =  +  a(U,  x)S  +  b(U,  x)  (3.12) 

with  boundary  conditions  (3.6).  The  functions  a  and  b  are  determined  in  terms  of  G*,  Hi 
and  H2.  We  multiply  (3.12)  by  pI5|P-25  ,p  >  2,  integrate  by  parts  over  [0, 1]  x  [0,  t] ,  where 
t  <T,  and  use  (3.9)  to  arrive  at 

^  \S{x,t)\Pdx  <  \So{x)\Pdx+K3P  J^{l  +  \S(x,T)\P)dxdr.  (3.13) 

Integrating  (3.13)  and  taking  the  1/p  power,  yields 

(  r  |5(x,<)|P(ix)'/P  <  [  /  \So{x)\Pdx  +  K3pTY^^e^^^] 

Jo  Jo 

therefore,  letting  p  — +  00,  we  obtain  (3.11).  | 
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4.  Convergence  to  Equilibria 

Let  (v(x,t),  cr(x,  t))  be  a  classical  solution  of  (1.1)-(1.5)  defined  on  [0, 1]  x  [0,  oo).  The 
existence  of  such  solutions  was  discussed  in  Section  3.  Here,  we  study  the  asymptotic 
behavior  of  as  f  — >  oo. 

The  first  lemma  indicates  that  5  =  cr  +  Ui  +  /x  converges  to  its  equilibrium  value. 
Lemma  4.1.  Under  ike  assumptions  of  Theorem  S.l, 


lim  5(x,  f)  =  0 , 

t— ►OO 

(4.1) 

uniformly  for  x  €  [0, 1]. 

Proof.  The  proof  is  a  consequence  of  the  a  priori  estimates 

/  f  Sfdxdr  <  C , 

Jo  Jo 

(4.2) 

r  r  S^dxdr  <  C  , 

Jo  Jo 

(4.3) 

f  5j(x,  t)dx  <  C ,  0  <  f  <  oo  , 

Jo 

(4.4) 

where  C  is  a  positive  constant  depending  only  on  data.  Once  these  have  been  established, 
the  identity 


/  S^{x,t)dx  =  f  (  S^{x,T)dxdT 
Jo  Jt-i  Jo 

+  2  f  f  f  S{x,r})St{x,Tj)dxdTjdT , 
Jt-l Jr  Jo 


(4.5) 


together  with  (4.2)  and  (4.3),  readily  yield 


lim 
<— *00 


(4.6) 


Then  (4.1)  follows  from  (4.4),  (4.6)  and  the  inequality 

S\x,t)  <  2{f  Sl{x,t)dxYI-‘  .  (4.7) 

It  remains  to  prove  (4.2)-(4.4).  By  (1-1),  (4.4)  is  contained  in  (2.7).  Also,  (4.3)  follows 
from  (2.8)  and  (2.7).  To  prove  (4.2),  observe  that  a  satisfies  the  equation 


<ru  +<7t  =  y'(Ur)Vrt 


(4.8) 
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whence,  one  derives  the  identity 


1  di 

2di  Jq  +  al{x,t)dx 

=  /  g'{vj:{x,t))vxt{x,t)(Tt{x,t)dx . 
Jo 


(4.9) 


Using  (2.7)  and  (2.17),  (4.9)  implies 


rf^‘ 


dxdt  <  C . 


By  (1.3), 


St  =  Wxf  +  at 


(4.10) 


(4.11) 


Combining  (4.11)  with  (4.10)  and  (2.7)  we  arrive  at  (4.2).  | 

Use  of  (4.1)  enables  us  to  identify  the  limiting  behavior  of  solutions  of  (2.11)  as  t  — >■  oo. 
The  following  result  is  analogous  to  Lemma  5.5  in  Pego  [14]. 

Lemma  4.2.  Let  s(x,«)  €  C^[0,oo)  he  the  solution  of  (2.11),  where  S(x,»)  is  continuous 
and  satisfies  (4.1),  0  <  x  <  1.  Then  s(x,»)  converges  to  s°°(x)  as  t  oo  and  s°°(x) 
satisfies 

s°°(x) g(s°°(x))  =  f  X  ,  0<x<l.  (4.12) 

Proof.  Let  x  6  [0, 1]  be  fixed.  Set  s(x)  =  liminft_oo  .s(x,  <)  and  s(x)  =  limsup(_op  s(x,t). 
By  virtue  of  (2.16),  — oo  <  s(x)  <  s(x)  <  oo. 

Assume  that  s(x)  <  s(x)  and  choose  Sq  €  (s(x),s(x)).  Then,  there  axe  sequences 
{t“},  tending  to  infinity  as  n  -+  oo,  such  that 


^((i,^)  ^  st(x,t~)  <  0,  s(x,t^)  =  So 

and,  by  (2.11), 

■Sf(3:,i^)  +  5o  +y(3o  -  S(x,tt))  =  fx. 
Let  n  — *■  CO  in  (4.14)  and  use  (4.13)  and  (4.1)  to  conclude 

So  +g(so)  =  fx. 


(4.13) 


(4.14) 


(4.15) 


Thus,  if  s(x)  <  s(x)  any  sq  €  (s(x),s(x))  satisfies  (4.15).  Since  the  solution  set  of  (4.15) 
is  not  connected,  this  is  impossible.  Hence  s(x)  =  s(x)  =:  s®°(x). 

Next,  we  claim  that  s®®(x)  =  lim{— oo  ^(x,  f)  satisfies  (4.12).  Indeed,  given  a  fixed 
X  in  [0, 1],  there  are  sequences  {in}  and  {r„}  ,f„  <  t„  <  tn  +  1  and  —*■  oo,  such  that 
(st(x,rn)|  =  |s(x,i„  +  1)  —  s(x,  t„)|  <  1/n.  Evaluating  (2.11)  at  (x,r„)  and  letting  n  — oo, 
we  deduce  that  5°°(x)  satisfies  (4.12)  for  any  0  <  x  <  1.  | 

Let  (v(x,t),  cr(x,t))  be  a  classical  solution  of  (1.1)-(1.5)  on  [0,1]  x  [0,  oo).  Recalling 
the  definition  of  s  in  (2.10),  Lemma  4.2  implies 


cr°°(x)  =  lim  <7(x,  t)  =  s°°(x)  —  /  x  . 

t— *oo 


(4.16) 
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Also,  combining  (1.3),  (4.1)  and  (4.16), 


v^(x)  ;=  limt;x(x,<)=  ]im  (S{x,t)  —  s(x,t))  =  ~s’^(x) 

t— ^oo  t— *oo 


(4.17) 


and 

5‘*(x)  =  uf  (x)  +  a°°(x)  +  /  X  =  0 . 
Finally,  noting  that  by  (1.4) 


Vx{x,t)dx , 


u°®(x)  is  Lipschitz  continuous  and  satisfies 


Yx)  :=  lim 

^  '  t— oo 


v(x 


(4.18) 


(4.19) 


(4.20) 


We  conclude  that  any  solution  of  (1.1)-(1.5)  converges  to  one  of  the  steady  states.  If 
0  <  /  <  m,  then  there  is  a  unique  smooth  steady  state  which  is  the  asymptotic  limit  of  any 
solution.  However,  if  m  <  /,  then  there  are  multiple  steady  states  and  thus  a  multitude 
of  possible  asymptotic  limits.  In  Sec.  5,  we  identify  stable  steady  states. 

Observe  that  in  case  w{^)  is  monotone  the  above  cirguments  yield  that  every  solution 
converges  to  the  unique  steady  state.  Moreover,  the  above  results  can  be  routinely  gen¬ 
eralized  to  the  case  that  the  function  w(^)  has  multiple  loops  but  the  graph  of  w  has  no 
horizontal  segments. 


5.  Stability  of  Steady  States 

The  scope  of  this  Section  is  to  study  the  stability  of  velocity  profiles  with  kinks.  To 
fix  ideas,  let  (u(x),a(x))  be  a  steady  state  of  (1.1)-(1.4)  such  that  tJ(x)  has  a  finite  number 
of  kinks  located  at  the  points  Xi, . . .  ,x„  in  (0, 1);  accordingly,  i;i(x)  eind  a(x)  have  a  finite 
number  of  jump  discontinuities  at  the  same  points.  Recall  that,  if  we  set  u(x)  =  — Uj;(x), 

Uj(tl(x))  =  ^X,X€[0,  l],X^Xi,...,Xn  (^•^) 


and  o(x)  =  g(— tl(x)). 

Given  any  smooth  initiad  data  (uo(a:),  o-o(x)),  there  is  a  unique  smooth  solution  (u(x,  <), 
cr(x,t))  of  (1.1)-(1.5).  As  t  oo,  (v(x,t),cr(x,t))  converges  to  one  of  the  steady  states, 
not  a-priori  identifiable.  We  now  restrict  attention  to  initial  data  (t;o(x) ,  cro(a:))  such  that 
the  values  of  (uoi(a^)  j  <^o(^))  ^ure  close  to  (ux(x),a(x))  except  on  the  union  U  of  small 
subintervals  centered  around  the  points  xj, . . .  ,x„.  U  can  be  thought  of  as  the  location  of 
transition  layers  separating  the  smooth  brauiches  of  the  steady  state.  Roughly  speaking, 
it  turns  out  that  the  steady  state  is  “stable”  under  smooth  perturbations  that  are  close  in 
energy,  provided  (t;(x),a(x))  takes  values  in  the  monotone  increasing  parts  of  w(^).  More 
precisely; 
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Theorem  5.1.  Let  (^v(x'),a(x))  be  a  steady  state  solution  as  described  above  and  satisfying 

u;'(u*(x))  >  Co  >  0,  X  6  [0,1],  x^xi,...,Xn  (5.2) 

for  some  positive  constant  cq.  If  the  measure  ofU  is  siifficiently  small,  there  is  a  positive 
constant  depending  onU  such  that,  if  6  <  Sq,  then  for  any  initial  data  (^0(2^),  o’o(x)) 
satisfying 

sup  |5o(x)|  <  6,  (5.3) 

0<z<l 


v'l{x,Q)dz  <'^6'^ 


(5.4) 


|uoi(a:)  -  Uj(x)|  <  ^  ,  x€[0,l]\W  (5.5) 

the  corresponding  solution  (t;(x,  t),  <t(i,<))  approaches  the  steady  state  (u(x),a(x))  as 
t  —*  00,  in  the  sense, 

Vx(x,t) -^Vj:(x),  (5.6) 

a(x,<) a(x),  (5.7) 

for  all  X  G  [0, 1]  \  l(. 

Proof.  To  simplify  the  exposition,  we  prove  the  theorem  for  the  case  that  u(ar)  =  —Vx{x) 
has  one  single  jump  discontinuity  located  at  xq,  m  <  /xq  <  M ,  cind  for  U  —  (xq  —e,  xq  +e) 
for  some  small  e.  Minor  modifications  are  needed  to  account  for  the  general  case. 

The  proof  is  based  on  exploiting  the  energy  identity  (2.6),  which,  upon  setting  u(x,t)  = 
—Vc(x,t)  yields  the  inequality 


~  J  v‘^(x,t)dx  +  J  [(W(u(x,t))  -  x/u(x,<))  -  $(x)]dx 

<7:  [  vf{x,Q)dx+  j  [(W(uo(a:))  -  2/uo(a^)  -  ^(aj)]<^a^- 
^  Jo  Jo 


(5.8) 


The  fimction  $(x)  is  associated  with  the  particular  choice  of  u(x),  it  identifies  a  basin  of 
attraction  of  u(x)  and  is  constructed  below. 

For  each  fixed  x,  the  function  u(x)  occupies  the  bottom  of  one  of  the  wells  of  the 
potential  I'V’(u)  —  xfu,  in  the  left  well  if  x  <  Xo  and  in  the  right  well  if  x  >  xq. 

Let 

e(x)  :=  W(u(x))  —  x/u(x)  (5.9) 

be  the  corresponding  value  of  the  potential  at  this  x,  eind  let 


E  :=  inf  inf  ( PF(0  -  x/0  >  —xf 


(5.10) 


be  a  lower  bound  of  the  potential  valid  for  all  x  G  [0,  Ij.  The  function  $(x)  is  now  defined 
by 

if|x-Xo|>£ 

if|x-xo|<6.  ^  ’ 
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Fig.  4:  Graphs  of  W(()  —  xf^,  e(x)  and  ^(x)  for  x  fixed,  x  <  xq  —e 
and  fx  above  the  Maxwell  line. 

Moreover,  set 


J  €  R  :  W’(^)  -  x/^  >  e(x)  for  — oo  <($<(},  x  <  xq 

\inf{(;eR:W{{)-xfl>e{x)  for  C  <  <^  <  oo},  x  >  xq 


and  observe  that  ^(x)  equals  +oo  or  — oo  in  case  u(x)  is  positioned  at  the  lower  of  the  two 
wells,  while  ^(x)  is  finite  otherwise.  The  functions  W’(^)  —  x/^,  e(x)  and  ^(x)  used  in  the 

construction  of  $(x)  axe  shown  in  Fig.  4  for  some  x  fixed,  with  x  <  Xq  —  e  0ind  fx  above 

the  McLxwell  line. 

The  following  lemma  identifies  the  region  of  attraction  of  the  steady  state  solution. 
Lemma  5.2.  There  is  a  positive  constant  po  such  that,  given  p  <  po,  there  are  positive  func- 
tzons  r and  r^^^p^,  vjzth  '^Tn^p'^  and  the  property^  For  any  r  ^  {p),rM{p)), 

if 

sup  |5o(x)|<p,  (5.13) 

0<r<l 

i^(0)  •=\  v^{x,0)dx  +  [^^(uoCx))  -  xfuQ{x)  -  ^x)]dx  <  ^p^  (5.14) 

and 

|uo(x)  -  u(x)|  <  r  -  2p,  xG[0,1]\W,  (5.15) 

then 

sup  |5(x,<)|  <  p  (5.16) 

0<I<1 

and 

|it(x,t)  -  u(x)l  <  r,  x6[0,l]\i/,  (5-17) 
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for  any  t  6  [0,  oo). 

Proof.  Let  t*  =  sup{to  €  :  supo<j.<i  15(x,t)|  <  p  for  0  <  t  <  <o}  be  the  first  time 

at  which  (5.16)  is  violated.  By  (5.13),  t*  >  0.  We  claim  that,  if  p  is  sufficiently  small, 
t*  =  +00.  Suppose  that  t*  <  +oo.  Then  for  any  to  <t*, 

sup  |S'(x,<)|  <  p,  0  <  <  <  <0  ,  (5.18) 

0<i<l 


while, 

sup  |5(x,t*)|=p.  (5.19) 

0<x<l 

We  now  refer  to  the  discussion  in  Section  2  on  invariant  intervals  of  uniform  length 
for  the  parametric  family  of  ode’s  (2.11).  Set  s(x)  =  “(2:)  =  ~Vx(x).  Let  pq^Rq  and,  for 
p  <  po,  let  ri(p)  and  r2(p)  =  Rq  —  ri(p)  be  as  in  Corollary  2.2;  also  recall  that  ri(p)  — >  0 
as  p  — »•  0.  Set  rm(p)  '■=  max{2p,ri(p)  +  p}  and  rjv/(p)  :=  mm{ilo,-Ro  —  n(p)  +  p}. 
Taking  po  even  smaller  if  needed,  it  follows  that  r„i(p)  <  rA/(p)  <  Rq  for  all  p  <  po.  Fix 


f  €  (t'm(p),rM(p))-  Then  equations  (1.6)  and  (2.10),  (5.13)  and  (5.15)  imply 

|so(x)  —  s(x)[  <  j — p,  x€[0,  l]\if.  (5.20) 

Moreover,  Corollary  2.2,  together  with  (1.3),  (2.10)  and  (5.18),  yields 

|u(x,<)  —  i£(x)|  <  |s(x,t)  —  s(x)|  +  |5(x,t)l  <  r  <  i2o  ,  x€[0, 1]\W.  (5.21) 

By  virtue  of  the  definition  of  Rq  and  (5.12),  Rq  <  ({x),  for  all  0  <  x  <  1 , x  7^  xq. 
Thus,  (5.11)  and  (5.21)  guarantee  that 

W(u(x,f))-x/u(x,t)  >  $(x),  0<x<l,0<t<<o.  (5.22) 

To  conclude  the  proof,  combine  (5.8)  with  (5.22)  to  obtain 

^  /  v1{x,t)dx  <  S(0),  0<t<to.  (5.23) 

2  Jo 

Then  (5.23)  and  (2.8)  together  imply  that 

sup  |5(x,fo)|  <  \/25(0)  (5.24) 

0<x<l 


for  any  to  <  <*.  But  (5.24)  and  (5.14)  together  contradict  (5.19).  Thus  t*  =  00  and  (5.16) 
holds  in  [0,oo).  Moreover,  (5.21)  yields  (5.17).  | 

A  remark  is  appropriate  at  this  point  to  guarantee  that  the  assumptions  (5.13)-(5.15) 
are  not  vacuous  but  follow  from  a  simpler  set  of  zissumptions  such  as  (5.3)-(5.5).  Recall 
that  the  constants  po  and  Rq  required  in  Lemma  5.2  depend  on  u(x),  and  the  form  of  w(^) 
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but  not  on  the  width  e  of  the  layer  in  the  initial  data.  If  5  <  min{p,T—2p},  then  (5.3)-(5.5) 
imply  (5.13)  and  (5.15).  They  will  also  imply  (5.14)  provided 

f\w{uo(x))  -  xf(uo(x))  -  ^x)]dx  <  jp"^ .  (5.25) 

Jo  ^ 

To  satisfy  (5.25)  two  things  axe  needed:  that  6  be  chosen  small  enough  and  that  e  not 
exceed  a  certain  value.  Moreover,  the  smaller  e  is  chosen,  the  larger  6  can  be  talcen. 

Lemma  5.2  guarantees  that  steady  states  (y^W)  satisfying  (5.2)  are  stable  under  pertur¬ 
bations  that  are  close  to  (v,a)  in  energy.  Actually  more  is  true.  It  follows  from  Lemma  4.2 
that  Vx(x,t)  —*  v^(x)  and  cr(x,t)  — >■  o’“(x)  for  a.e.  x  6  [0, 1]  eis  t  — ♦  oo.  Using  (5.17), 

lvx(x,t) -Vx(x)l  =  lu(x,t)  -  ic(x)l  <  r,  x6[0, 1]\7/,  (5.26) 

and  the  fact  that  the  only  solution  of  (5.1)  in  [t<(x)  —  r,u(x)  +  r],  for  r  <  Rq,  is  u{x)  itself, 
we  conclude  that  v^(x)  =  Vx(x)  and  o'°°(x)  =  '^(x)  for  x  E  [0, 1]  \  ZY-  | 


Appendix 

In  this  appendix,  we  use  the  results  of  Sec.  3  to  obtain  global  existence  results  for 
equations  governing  motions  of  incompressible  viscoelastic  fluids  in  simple  planar  shear. 
We  consider  a  class  of  one-dimensional  models  popular  in  polymer  rheology  based  on 
differential  constitutive  laws  which  satisfy  the  principle  of  frame  indifference;  we  refer 
to  [16]  for  a  general  formulation. 

1.  Johnson-Segalman  /  Oldroyd  Models  For  a  class  of  models  with  m  -  relaxation 
times,  derived  from  three-dimensional  differential  constitutive  relations  due  to  Johnson 
&c  Segalman  [6]  (with  exponential  memory  functions)  and  Oldroyd  [12],  the  equations  of 
motion  take  the  following  form  in  one  space  dimension  (see  [10]  for  a  detailed  formulation): 

pvt  =  Sx  , 
m 

<Tjt  +  Xj(Tj  =  {Zj  +  pj)  Vx  , 

Zjt  +  XjZj  =  —  (l  —  Gj)  (TjVx  ,  J  =  1)  —,Tn  , 

where  p  is  the  fluid  density,  v  is  the  fluid  velocity,  aj  represents  a  contribution  to  the  toteJ 
shear  stress  due  to  the  polymer,  rjVx  is  the  Newtonian  contribution  to  the  total  shear  stress, 
and  Zj  represents  a  contribution  to  the  principal  normal  stress  difference;  rj,  pj,aj 
are  positive  constants  where  —1  <  aj  <  1;  Xj  are  relaxation  rates,  pj  are  elastic  shear 
moduli,  and  Cj  are  slip  parameters. 

The  initial-boundary  value  problem  governing  Poiseuille  shear  flow  between  parallel 
plates  at  z  =  ±1  with  the  flow  symmetric  about  the  centerline,  and  driven  by  a  constant 


-20- 


pressure  gradienc  consists  of  the  system  (A.l)  on  [0,1]  x  [0,  oo),  coupled  with  boundary 
conditions  (1.4),  and  initial  conditions 

u(x,  0)  =  Vo(x),  <Tj(x,  0)  =  o-jo(x),  Zj(x,  0)  =  Zjo(x) ,  {A.2) 

on  0  <  X  <  1,  j  =  l,...,m;  it  is  assumed  that  the  compatibility  conditions  uo(l)  = 
0,  t;q(0)  =  0,  and  crjo(O)  =  0,  j  =  l,...,m,  are  satisfied.  In  what  follows,  we  assume 
that  the  initial  data  satisfy  the  hypothesis  of  Theorem  3.5.  We  note  that  when  r/  =  0, 
it  is  readily  shown  that  (A.l)  is  hyperbolic  if  speeds 

±{i  Y^Y-i(Zj  +  and  0  (repeated  m  +  1  times). 

In  the  case  of  a  single  relaxation  time,  we  omit  the  subscript  j,  and  we  non-dimensiona- 
lize  the  variables  by  scaling  distance  by  /i,  time  by  A“^,  and  stress  by  fi.  Furthermore,  if 
we  replace  a,  u,  and  /  by  d  ;=  (1  —  u  ;=  (1  —  and  /  :=  (1  —  a^)^/^/, 

respectively,  then  the  parameter  a  disappears  from  Eqs.  (A.l).  Since  no  confusion  will  arise, 
we  omit  the  caret.  There  are  two  essential  dimensionless  parameters:  o:  :=  ph'^X^/n,  a  ratio 
of  Reynolds  number  to  Deborah  number,  and  e  :=  t]\/ fx,  a  ratio  of  viscosities.  The  special 
case  of  (A.l)  with  j  =  I  leads  to  the  initial-boundaxy-value  problem  consisting  of  the 
system  (1.10)  on  [0, 1]  x  [0,  oo),  with  boundary  conditions  (1.4)  and  initial  conditions  (A.2), 
j  =  1. 

To  establish  a  global  existence  result  for  the  system  (A.l),  with  boundary  condi¬ 
tions  (1.4)  and  initial  conditions  (A.2),  we  observe  that  (A.l)  is  endowed  with  the  identities 


=  -2A^l(l-a>?+(Zj  + 


\2  _ 

2  ^  4 


(A.3) 


for  j  =  1,  ...,m,  and  any  x  €  [0, 1].  For  each  j,  the  right  side  of  (A.3)  is  strictly  negative 
in  the  exterior  of  the  ellipse 


r,:=  {{aj,Zj):{l-ajyj+iZj  + 


N_\2  _  ^7  •] 

2  ^  4  ^ 


(A.4) 


Consider  the  closed  sets 

Q,j  :=  {(o-j,  Zj)  :  (1  -  apcr|  +  (Zj  -h  pjY  <  Cf}  , 


(A.5) 


j  =  1,  ...,Tn,  where  Cj  >  0  are  constants.  H  Cj  >  =  1, ...,  m,  the  ellipses  bounded  by 

Fj  are  (properly)  contained  in  the  sets  Q,j  for  each  j.  Letting  U  :=  (cri,  Zi,  ^m) 

defining  :=  IJ^i  above  construction  implies  that  if  Cj  >  iJ,j{j  =  l,...,m),  U  € 

L=°([0, 1]  X  [0,T])  for  every  T  >  0,  where  the  a  priori  bound  depends  only  on  the  initial 
data  and  not  on  T.  Thus  Theorem  3.5  yields  that  solutions  of  the  initial-boundary  value 
problem  (A.l),  (1.4),  (A.2)  exist  globally  in  time. 

Finally,  we  remark  that  the  same  approach  can  be  applied  to  obtain  global-time 
existence  for  shezuring  flows  of  non-Newtonian  fluids  governed  by  differential  constitutive 
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relations  due  to  Phan-Thien  zuid  Tanner  [15]  .  In  the  case  of  one  relaocation  time,  the 
relevant  system  in  dimensionless  form  (compare  with  (1.10))  is 

OCVt  —  f  i 

at  +  =  (Z  +  l)vx  ,  (A.6) 

Zt  +  e^^Z  =  -crvx  , 

where  a,  e,  u,  a  have  the  same  meaning  as  in  (1.10)  and  where  the  constant  <5  >  0.  To  estab¬ 
lish  a  global  existence  result  for  the  initial-boundary  value  problem  associated  with  (A.6), 
we  use  the  same  approach  as  for  (A.l).  For  the  application  of  Theorem  3.5,  the  identity 
for  (A.6),  corresponding  to  the  identity  (A.3)  with  j  =  1,  now  takes  the  form 

+  (Z  +  1)^1  =  +  (Z+  -  i] ,  (.4,7) 
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